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Abstract 

We consider Grover's search algorithm on a model quantum computer implemented 
on a chain of four or five nuclear spins with first and second neighbor Ising interactions. 
Noise is introduced into the system in terms of random fluctuations of the external 
fields. By averaging over many repetitions of the algorithm, the output state becomes 
effectively a mixed state. We study its overlap with the nominal output state of the 
algorithm, which is called fidelity. We find either an exponential or a Gaussian decay 
for the fidelity as a function of the strength of the noise, depending on the type of noise 
(static or random) and whether error supression is applied (the 27r/c-method) or not. 

PACS: 03.67.-a,3.67.Lx, 03.65.Ta, 03.67.Dd, 03.67.Hk 

1 Introduction 

It is well known that Grover's search algorithm for searching an item in a quantum data base 
of length N = 2"^, where n is the length of the quantum register made up of n-qubits, provides 
a quadratic speed up with respect any classical search algorithm [H [2]. Some applications 
of this algorithm have already been done [3], HI O [6], [71 [8] , and its physical implementation 
has been made for few qubits P, [lOl [TT], [12] . Numerical simulations of this algorithm have 
also been made for a 4-qubits solid state quantum computer which is made up of a chain of 
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four Ising paramagnetic nucleus of spin one-half interacting with the first and second spins 
neighbors [13]. This simulations use RF-pulses to implement quantum gates [2], and the 2iTk- 
method [15] to control the non resonant transitions in the system. At the moment, there has 
not been a physical realization of this quantum computer. However, since the full Hamiltonian 
can be explicitly written for this system, several useful analytical and numerical studies have 
been made [TBI [T7] where the results can be extrapolated to other solid state spins quantum 
computers. On the other hand, one important aspect of the simulation or implementation of a 
quantum algorithm is the effect that the noise can cause to the performance of this algorithm 
which may arise from error in the quantum parameters or from the interaction with the 
environment, or both. This noise may disrupt the algorithm (decoherence) and, eventually, 
destroy it. Previous studies of the effect of noise on a realisation of Grover's quantum search 
algorithm have been based on description in terms of abstract gates, which have been turned 
into noisy gates [18], [191 ISOl [211 [22] and/or the introduction of dissipation and decoherence 
within the Born-Markov approximation [23]. By contrast, the present model is based on a 
definite quantum system, with a specific Hamiltonian, where unitary errors and noise have a 
clear physical interpretation. In our model, unitary errors arise from non-resonant transitions, 
whereas noise arises from fluctuations in the external fields - either the constant magnetic 
field or the radio-frequency pulses. 

In a previous paper [13], we studied unitary errors in the realisation of Grover's quantum 
search algorithm on our quantum computer model. The present work is devoted to the study 
of noise in essentially the same setting. In Sec. [21 we introduce the physical model quantum 
computer. In Sec. [31 we explain how the quantum search algorithm is implemented with the 
help of certain radio-frequency (RF-) pulses. Sec. [H describes the implementation of noise, 
by letting the Rabi frequency (Rabi-frequency noise), or the Larmor frequencies (Larmor- 
frequency noise) fluctuate. The effect of noise on the quantum computation is quantified 
using the fidelity, first introduced in [2l] in an attempt to transfer the concept of exponential 
instability from classical to quantum mechanics. For a recent review on fidelity, see |25j . 
Sec. [SI presents our actual numerical results, and we summarize our findings in Sec. [B 



2 Hamiltonian of the model quantum computer 



The model of the system can be found in reference [12], and the explicit Hamiltonian for 
simulation of the Grover' search algorithm is given by [T7] 

H = Ho + W{t) , (1) 
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where one has that 

^32 1 

^H, = -J2^,F,-2JJ2 n 4\i - 2 J' Yl 4^2 , (2) 

k=0 k=0 k=0 

and 

3 

W{t) = - ^ E {e'^-'^-^n + e-'(-*+-)4-) , (3) 

fc=0 

where Uk is the Larmor frequency of the /c'th-spin, J and J' are the couphng constant to first 
and second neighbor spins, Q is the Rabi frequency, u and ip are the frequency and phase of 
the radio-frequency pulse. The operators and are single qubit operators acting on 

the kth-qubit of the registers {|q;3Q;2Q;iq;o)} with = 0, 1 for z = 0, 1, 2, 3. These 2^-states are 
eigenstates of the Hamiltonian Hq, and the action of the operator on is associated qubit is 
as follows 

= ^^l"'^) • (4) 

The operator = (I^Y = |Ofc) (1^0 represents the descend operator for the qubit \ak), and 
the energy difference between two eigenstates of Hq is 

= h{wk + f^J + i^J') , /i,!^ G {-2,-1,0,1,2} . (5) 

Any quantum gate is realized within the interaction picture, applying defined RF-pulses with 
rectangular shape, choosing the frequency lv in resonance with the appropriated transitions, 
and given the time (r) duration of the pulse. The corresponding evolution operator for the 
whole quantum system, is denoted by R'^'^{ip,9), where 9 = Qt, and the RF- frequency is set 
equal to w = + fiJ + uJ' (resonant situation). 

Degeneracies in Ho At a first glance, one would expect accidental degeneracies in the 
Hamiltonian Hq to cause problems in the realisation of any quantum algorithm, since a 
given microwave field could drive unwanted transitions involving degenerated states. Such 
degeneracies can be avoided requireing the Larmor frequencies Wk to increase in powers of 
two: Wk oc 2^. Without interactions J = J' = 0, this produces an equidistant spectrum; 
for J' <^ J <^ this guaranties a spectrum without any degeneracies. Such a choice of 
the Larmor frequencies would destroy the scalability of the present system as a model for 
quantum computation. 

^Note, in standard spin algebra, we would expect /+ = + i/^ = |1)(0| being the ascend operator. 
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However, as it turns out degeneracies are not really a problem. This is due to the fact 
that the interaction operator couples states la) and only if they differ in no more than 
one qubit. In the interaction picture, W{t) becomes: 

Wo,p{t) = -^\ e"*(^^-*+^) : a - /5 = (0, . . . , Ifc, . . . , 0) , Ak = w - {wk + ^iJ + uJ') , 
[o : else 

(6) 

where the notation a — (3 = (0, . . . , 1^, . . . , 0) means that the multi-indices a, (3 are equal 
except at one position (qubit k) , where the qubit in the a-state is excited and the qubit in the 
/3-state is not. For an accidental degeneracy caused by the unfortunate choice of the Larmor 
frequencies, there must exist transitions such that + fiJ + uJ' = Wi + fi'J + u'J' for some 
k ^ I. But this can be savely avoided by choosing the spacing between neighboring Larmor 
frequencies much larger than the coupling constants: \wk — Wk-i\ <^ J <^ J' . Since this 
allows for a linear increase in the Larmor frequency as a function of the number of qubits, 
the scalability of the quantum computer model is not affected. 



Near resonant approximation Assume that the frequency w of the RF-pulse matches 
a particular transition so well that all other transitions can be neglected. For a rectangular 
RF-pulse as described above, with detuning A. = w — Wk — /U J — uJ\ the evolution operator 
in the basis of the two states involved read: 

_ X /cos - ^ sin %^ f e^^ sin ^ \ 

V e--^^^") \ f e-^ sin ^ cos + sin i ' ^ ' 



where VL^. = + A^. For A = 0, we obtain the resonant transitions: 

/ cos^/2 ie*^sin^/2' 
Vie ^ sm6'/2 cos 6/ 2 

3 Implementation of Grover's search algorithm 

For our case, Grover' search algorithm [26] requieres three qubits to prepare the "enquiry" 
states and a single qubit (" ancilla" ) required by the oracle to communicate its answer. Starting 
from the ground state, IO3O2O1O0), a superposition state with the qubits = 0, 2, 3 is generated 
by applying the corresponding Hadamard gates, 

H^'^ = HSH2H0 , (9) 
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to this state, where is the Hadamard gate acting to the kth-qubit. Then, one apphes the 
Grover operator, 

G = OaH^^^ SoH^'"^^ , (10) 

twice, since for = 8 = 2'^ registers, the probabihty to find the searched state is maximum 
for [vr a/8/4 — 1/2] ^ 2 apphcations of Grover operator (steps). The operator 5*0 represents 
the phase inversion of the searched state, and the operator Oa represents the oracle where a 
is the index (state) to be found in decimal notation. The implementation of these operatos 
by pulses is found in [13], 

^0 = {u,,.=-^,,Kih^)} {n,,=_i,i^r(f,f)} {n,,=-M^r(^>^)} • (n) 

H2 = {n,,=-2,o,2^r(f,vr)} {n,,=-w^r(f,f)} {n,,=-M^r(^,^)} • (12) 
H3 = {n„.=-i,i^r(f,^)} {n„.=-i,i^r(f,f)} {u,,.=-,,^Rri^,^)} ■ m 

So = RI'~\0, 27r) R['\0, 27r) i??'~^(0, 2tt) Rf'~\0, 27r) R^^'^O, 27r) , (14) 

Oo = i??''(0,27r) Os = R'{-\0,27r) (15) 

O5 = R^^'\0, 27r) On = i?r''"'(0, 27r) . (16) 

One applies 70 pulses to represent the operator H^^^ , for the superposition state, and 146 
pulses to represent Grover operator to have a total of 362 pulses to implement Grover' search 
algorithm. 

The duration of the RF-pulses used depends on the rotation angle 6 and on the Rabi 
frequency: Tpuisc = 0/fl. Further below, we will measure time in units of tph = 7i/{2Q), the 
duration of a 7r/2-pulse. 

4 Implementation of noise 

The Larmor frequencies Uk for /c = 0, 1,2,3 depend on a static magnetic field in z-direction 
with a strong field gradient: Uk = 7-Bfc, where is the strength of the static magnetic 
field at the location of the fcth-paramagnetic nucleus of spin one-half. The Rabi-frequency 
Q is determined by the strength of the radio-frequency (RF-) pulses which is assumed to be 
constant during a pulse: Q = 'yBt, where 7 is the gyromagnetic ratio and Bt is the amplitude 
of the transversal RF-pulse. The environment or imperfections in the magnetic field produce 
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noise which, in turn, will produce variations in the above parameters. To study the effect of 
these variations, we consider random detunings of the form 



i^k = '^ko + eLi (17) 

and 

n = n^ + eR^, (18) 

where ujko for /c = 0, 1, 2, 3 and flo are the Larmor and Rabi frequencies without noise, ^ is a 
random Gaussian variable, centered at zero, with unit variance. The parameters ei and en 
determine the amplitude of the noise. 

The two types of noise produce different types of errors as can be seen from the evolution 
operators in the resonant approximation, Eqs. ([7]) and ([8]). In the case of noise on the Rabi 
frequency, the evolution operator remains of the form ([8]), only the angle 6 suffers a change: 

6 ^6' = 9 + An, An = eR^r, (19) 

where r is the duration of the RF-pulse. In the case of noise on the Larmor frequencies, the 
evolution operator gets distorted in a more complicated way: With 6 = Qt and = ei^r 
we find: 

/eiA./2 \ /cose/2 -i^ sin 0/2 ie'^ sine/2 

u ^ / y le ^ sm6'/2 cos6'/2H — ^ sm6'/2 

'(20) 

up to second order in = ei^r. In both cases, the strength of noise must be compared to 
the smallest energy (frequency) scale in the system. This is J' the coupling between second 
nearest neighbors. 

Besides the amplitude also the frequencies of the noise fluctuations are important. It 
makes a difference whether the variation of the Larmor frequencies (Rabi frequency) happen 
on a time scale which is of the order of the duration of the whole algorithm (low frequency 
noise) or of the duration of a single RF pulse. To investigate this issue, we consider two 
types of noise (errors): Static noise, where the respective frequencies are detuned for each 
realization, but kept fixed during the whole algorithm. Random noise, where the respective 
frequencies are detuned anew for each RF-pulse. The latter simulating noise which fluctuates 
on a time scale of the order of the inverse Rabi frequency. 



Fidelity In order to measure the effect of noise on the realisation of Grover's quantum search 
algorithm, we use the fidelity (or quantum Loschmidt echo) [2H [25] . We perform averages 
over TT-rep repetitions (where rirep is of the order of 100) of the algorithm. This effectively turns 
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the final state of the quantum computer into a mixed state ftcai- The fidehty is computed by 
projection of ftoai on the outcome of the ideal realisation of the algorithm: 

1 " 

i^end = (^ideall ftcal l^ideal) , ftcal = I ^leil) (^ieil I ) (21) 

^ r=l 

where |^I/fcL) final states obtained from the algorithm taking into account the random 

detuning of the respective frequencies. 

When analysing the final fidelity, there are two options: we can compare the final state 
after the evolution in the presence of noise and unwanted transitions with the final state 
obtained from an error-free execution of the quantum protocol, or we can compare the final 
state with the desired target state. These two options are not the same, because the Grover 
protocol arrives at the desired target state only with a final probability, which is close to but 
still different from one [27] . Both options yield very similar results for the behavior of fidelity. 
For the results presented in this paper, we chose the first option. 

5 Numerical method and results 

We chose J = 10 (one order of magnitude smaller than the spacing between the Larmor 
frequencies) and J' = 0.4 (one order of magnitude smaller than J). In this way, made 
sure that the dominant source for errors were the near resonant transitions with detunings 
of the order of J'. These parameters are all frequencies - see Eq. ([2]) - given in units of 
27T X MHz. Note however that since the model depends only on their relative values, we 
may also choose any other unit. Unless stated otherwise, we choose the Rabi frequency as 
Qo = 0.1008 according to the 27rA;-condition, such that near resonant transitions due to the 
second-neighbor coupling are efficiently suppressed. 

The evolution of the quantum state in the qubit-register is calculated numerically by 
solving the time-dependent Schrodinger equation in the interaction picture. For that we have 
to solve a first order ordinary differential equation, where the state vector is represented by 
N = 2^ complex amplitudes. In view of the rapidly oscillating right hand side (the interaction 
matrix), any sophisticated adaptive stepsize routine is of limitted use. We therefore applied 
a simple 4'th order Runge-Kutta algorithm to solve the differential equation, where we fixed 
the stepsize At according to the largest Larmor frequency Wn which is the largest frequency 
in the problem: At = 7r/(2ty„). 

The computation time depends mainly on the time needed to evaluate the right hand 
side of the differential equation. We have to calculate the components of an A^- dimensional 
vector, and for each component we have to sum over the n non-zero matrix elements of the 



7 



interaction matrix. This means that the time needed for the evaluation of the right hand side 
scales as nN = 2""*"^ with the number of qubits. In addition, more qubits require a larger 
pulse sequence to apply the Grover algorithm. Ultimately, the length of the pulse sequence 
would scale with \/N as explained for instance in Ref. [26]. This would mean that overall, 
the computation time scales as 2^"/^ . The simulation presented in this work are done on 
standard single processor machines. With our hardware a single run of the Grover algorithm 
with 4 or 5 qubits takes approximately 4 minutes computing time. 

Choice of the Larmor frequencies On the basis of the discussion of accidental degen- 
eracies in section [21 we expect that the Larmor frequencies may be chosen to increase linearly 
with the number of qubits. Thus, we normally set 

Wo = 50 , wi = 200 , W2 = 350 , = 500 (22) 

for the calculations with four qubits and in addition = 650 for the 5-qubit calculations. 
However, before settling for these values we compared the errors in different schemes (expo- 
nential increase, avoiding any degeneracies in Hq, quadratic increase, and the present linear 
increase). We found consistently, that the different schemes only lead to apreciable errors if 
at least one of the differences Wk — wi gets close to the energy scale of the coupling constants. 
Our choice in Eq. fl22l) is perfectly save in view of the values chosen for J and J'. 

5.1 Larmor frequency noise 

To start our study, we analyse single runs without any averaging. We compare calculations 
with and without noise, as well as exact calculations and calculations which apply the near 
resonant approximation. We use a quantum register of 3+1 qubits to simulate Grover's search 
algorithm. The qubits are arranged according to the string xx.x, where x represents a data 
qubit, while the . represents the auxiliary qubit which is used to implement the oracle and 
the conditional reflection. The algorithm consists of an initial sequence of Hadamard gates 
H^^^ followed by two Grover steps. 

In the simulation for Fig. [T]the algorithm searches for the state (string 00. 0). On the 
graph we indicate the time intervals where the quantum oracle is implemented with red bars. 
These intervals therefore mark the start of a Grover step. Furthermore, we mark the time 
intervals where the conditional reflection is implemented with green bars. The red solid line 
shows the fidelity as a function of time for the exact calculation without noise. During the 
algorithm, the fidelity drops with more or less constant rate from -F(O) = 1 to Fend ~ 0.92. 
This reminds us that the implementation of the quantum gates via RF-pulses as such already 
introduces errors. On the graph, time is indicated in units of thp the duration of a 7r/2-pulse 
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Figure 1: Fidelity as a function of time for single runs of Grover's quantum search algorithm. 
Numerical simulations: without noise and exact (red line), using the near resonant approx- 
imation (green line); with static Larmor frequency noise [e^ = 0.02) and exact (blue line), 
using the near resonant approximation (pink line). The red and green vertical bars indicate 
the position of pulses implementing the oracle (red bars) and the conditional reflection (green 
bars). Time is indicated in units of t^p the duration of a 7r/2-pulse. 

(see end of Sec. [3]). The green solid line shows the same calculation, but with the near resonant 
approximation being applied. It follows the exact calculation very well. The near resonant 
approximation neglects the far resonant, where the frequency missmatch is of the order of the 
Larmor frequencies. The fidelity is generally a bit larger than for the exact calculation. The 
agreement between the exact calculation and the near resonant approximation is completely 
lost as soon as the Larmor frequency noise (here it is static noise may be better described 
as random detunings) is switched on (blue and pink solid lines). The blue line shows the 
exact calculation for ei = 0.02 which only initially follows the curve for the fidelity decay 
without noise. More or less in the middle of the algorithm, the fidelity curve drops rather 
abruptly to a level of F{t) ~ 0.7, where it saturates. The near resonant approximation doesn't 
show this behavior. Instead it stays close to the curve which shows the fidelity decay without 
noise. We may conclude that the noise introduces additional errors primarily via such indirect 
transitions which are neglected in the near resonant approximation. In the simulations shown, 
there is no indication that the oracle or the conditional reflection play any special role in the 
loss of fidelity. 

Figure [2] shows the average fidelity as described by Eq. (!2T]) for static and random noise on 
the Larmor frequencies. The simulations are averaged over an ensemble of rirep independent 
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50 100 150 200 250 300 350 400 

Figure 2: The average fidelity as a function of time in the case of static (red hues) and random 
(blue lines) Larmor frequency noise. Exact numerical simulations and the near resonant 
approximations are shown in the same color. The latter largely over-estimate the fidelity and 
fail to reproduce its true behaviour. Red and green vertical bars again indicate those time 
intervals, where the oracle and the conditional reflections are applied. 

repetitions of the algorithm. In the case of the exact numerical simulations, n^ep = 50 (static 
noise) and 25 (random noise). In the case of the near resonant approximations, rifep = 1000 
in both cases. This figure serves to demonstrate that the behavior of the fidelity may be very 
different for different repetitions - at least for static noise (compare with Fig. [1]). We expect 
smaller variations within individual repetitions in the case of random noise, since there, the 
fidelity depends on a large number of random variables which tend to wash out any kind of 
extreme behavior. Interestingly, the simulations reveal that for static noise, fidelity is lost 
to a large extent during the oracle and the conditional reflection pulses. For random noise, 
the fidelity drop is more uniform. Finally, Fig. [2] again demonstrates the failure of the near 
resonant approximation in the presence of Larmor frequency noise. 

Static vs. random Larmor frequency noise 

Here, we analyse the average fidelity for the whole algorithm as a function of the noise 
amplitude. Averages are now taken over ensembles which contain typically about Uj-cp = 100 
repetitions. The statistical error (due to the finite size of the ensemble) is estimated as 
follows: We divide the Up fidelity values for each repetition into p sub-ensembles of equal 
size. This yields average fidelities Fi, i = 1, . . . ,p. We take care that these sub-ensembles 
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Figure 3: The fidelity Fend ^oi static (a) and random (b) noise on the Larmor frequencies. 
Exact numerical simulations: 4 qubits, target state (red error bars); 4 qubits, target state 
8 (green error bars); 5 qubits, target state 27 (blue error bars). The solid lines of the corre- 
sponding colour show a phenomenological fit with Eq. ^IM to the numerical data. 

are sufficiently large such that the group-averaged fidelities are Gaussian distributed. That 
allows us to estimate the variance of the full average: F = Yl^=i calculating the 

best estimate for the variance of the group averages: 

1 1 ^ 

var(F) = - var(F,) , var(Fi) = V(F, - F)^ . (23) 

p p — 1 ^ — ^ 

i=l 

We then use the standard deviation (j{F) = A/var(F) as an estimate for our statistical error. 

Fig. [3] shows the final average fidelity Fend of the whole algorithm as a function of the 
amplitude of the Larmor frequency noise. We compare the simulations with static noise 
(panel (a)) and random noise (panel (b)). In either case, we simulate the algorithm for four 
qubits and target states and 8, and also for five qubits and the target state 27. For a given 
noise model (static or random) the fidelity Fend behaves very similar in all cases. In panel (a) 
we observe a roughly exponential decrease for increasing ei- At large values of the fidelity 
saturates at approximately /bas ~ 0.175. In panel (b), random noise, we observe a somewhat 
faster decay. Again the fidelity saturates for large e^, but now the fidelity saturates earlier 
than in the case of static noise. Also, the saturation level is somewhat smaller: /bas ~ 0.08. To 
compare the behavior of the final fidelity Fond more quantitatively, we consider the following 
phenomenological model function: 

/(e) = /bas + (1 - /bas) exp [-(e/e2)' - e/ei] , (24) 
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with the fit parameters /bas? ^2; and ei. In this parametrization, €2 indicates the point where 
Gaussian decay becomes important. Similarly, ei indicates the point where exponential decay 
becomes important. For the three cases shown in Fig. [3](a), static noise, we obtain the 
following results: 





/bas 


£2 


ei 


4 qubits, 0(0) 
4 qubits, 0(8) 
5 qubits, 0(27) 


0.161(49) 
0.186(19) 
0.178(62) 


0.046(20) 
0.0359(40) 
0.11(22) 


0.0315(73) 
0.0510(80) 
0.0235(23) 



(25) 



where the statistical uncertainties of the fit values are given in parenthesis (they refer to the 
last two digits of the respective fit value). For the three cases shown in Fig. [3](b), dynamic 
noise, we obtain: 





/bas 


£2 


ei 


4 qubits, 0(0) 
4 qubits, 0(8) 
5 qubits, 0(27) 


0.104(12) 

0.0904(45) 

0.0555(83) 


0.0224(20) 

0.01870(55) 

0.0209(14) 


0.0303(44) 

0.083(16) 

0.0291(34) 



(26) 



Comparing the reduction of the fidelity for static and random noise, we find that random noise 
has a more destructive effect. This holds even at very small noise levels. This result may be 
surprising in view of opposite results in Refs. [T21[2H]- However, note that both articles focus 
on larger quantum registers with more qubits and correspondingly more Grover steps. This 
is possible because they work with abstract gates and gate errors. 

5.2 Rabi- frequency noise 

In this section we study the effects of amplitude noise on the RF-field. In the present approach 
this is equivalent to adding noise to the Rabi-frequency. As in the previous case, we consider 
static and dynamic noise. In the case of static noise, the fidelity depends on only one random 
variable, the detuning of the Rabi-frequency: f2 — > f2o + ^r^, where ^ is a random Gaussian 
variable with zero mean and unit variance. This detuning of the Rabi frequency changes at 
random for each repetition, while being kept fixed during the execution of the algorithm. In 
the case of dynamic noise, the detuning changes at random for each microwave pulse. As in 
the case of Larmor frequency noise, we will first study the behavior of fidelity as a function 
of time and only afterwards the final fidelity (for the execution of the whole algorithm) as a 
function of the noise amplitude. We will find that the near resonant approximation is accurate 
in the presence of Rabi frequency noise. In the study of the final fidelity -Fend? "we therefore 
use that approximation. 
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Figure 4: Fidelity as a function of time for single runs in the presence/absence of static Rabi 
frequency noise. Numerical simulations: without noise and exact (red line); using the near 
resonant approximation (green line); with noise, ejiC, = 0.005, and exact (blue line); using the 
near resonant approximation (pink line). Red and green vertical bars again indicate those 
time intervals, where the oracle and the conditional reflections are applied. 

Fig. m shows the fidelity as a function of time for single runs of Grover's search algorithm. 
We compare simulations with and without static Rabi frequency noise in analogy to Fig. [H 
where we showed a similar figure for the case of static Larmor frequency noise. Here, a single 
run with static Rabi frequency noise, rather means that the simulation is performed with a 
certain detuning: enC, = 0.005. Comparing these results with those in Fig. [H we find that 
a certain degree of fidelity loss requires a considerably larger amount of noise in the case 
of Larmor frequency noise (e^ = 0.02). Furthermore, we observe that the near resonant 
approximation is very accurate in the present case of Rabi frequency detuning. 

Fig. [5] shows the average fidelity for static (red line with error bars: exact numerical cal- 
culation; green line with error bars: near resonant approximation) and dynamic (blue line 
with error bars: exact numerical calculation; pink line with error bars: near resonant approx- 
imation) Rabi frequency noise. The error bars represent the estimate of the statistical error 
as explained in Eq. fl23|) . The exact calculations are performed with rirep = 100 repetitions, 
the near resonant approximations with n^cp = 1000. Note that the statistical uncertainty is 
quite large for static noise whereas it is very small for dynamic noise. Comparing to Fig. [2] 
we observe much less pronounced differences between oracle and conditional reflection on the 
one hand and the Hadammard gates on the other. In both cases, static and random noise, the 
difference between exact simulation and near resonant approximation is within the statistical 
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Figure 5: Average fidelity as a function of time for static and dynamic Rabi frequency noise 
with amplitude = 0.005. Numerical simulations: Static noise, exact (red line); near 
resonant approximation (green line). Dynamic noise, exact (blue line); near resonant approx- 
imation (pink line). 



Static vs. Random Rabi frequency noise 

From the results for the fidelity as a function of time, we saw that the near resonant approx- 
imation works very well in the case of Rabi frequency noise. Since it provides an enormous 
speed-up of the simulations, we will use it almost exclusively in the present section. 

Figure [6] shows the average fidelity Fend for the complete algorithm as a function of the 
amplitude of the Rabi frequency noise. In panel (a), we consider static noise, in panel 
(b), dynamic noise. We find a striking difference in the functional dependence of the fidelity 
between the two cases. For dynamic noise (panel (b)) we expect the fit function Eq 
to work well. However in the case of static noise (panel (a)) we clearly need a different fit 
function to account for the apparently algebraic decay. In an attempt to maintain the general 
structure of the fit function as much as possible, we choose the new fit function of the form 



The following two tables contain the results of the fits for static and dynamic Rabi frequency 
noise. 



error. 




(27) 
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Figure 6: The fidelity F^^d for static (a) and dynamic (b) Rabi frequency noise of strengtli 
cr. Numerical simulations using the near resonant approximation: 4 qubits, target state 
(red error bars); 4 qubits, target state 8 (green error bars); 5 qubits, target state 27 (blue 
error bars). The solid lines of the corresponding colour show the phenomenological fit with 
Eq. (1271) (panel (a)) and Eq. (124|) (panel b) to the numerical data. The black error bar near 
the red curve shows one case, where the fidelity has been calculated exactly (static noise, 4 
qubits, target state 0). 

(a) Static noise: For the fit function g{e) from Eq. (127|) one obtains the following fit values 
for the three cases considered: 





ghas 


£2 


ei 


4 qubits, 0(0) 
4 qubits, 0(8) 
5 qubits, 0(27) 


0.052(10) 

0.0230(98) 

0.0798(69) 


0.00429(23) 
0.00633(24) 
0.00289(14) 


0.0140(45) 

0.032(10) 

0.0064(14) 



(28) 



(b) Dynamic noise: Here, we use the same fit function /(e) from Eq. ( !24l) as in the case of 
Larmor frequency noise, and we obtained the following fit values: 





/bas 


£2 


ei 


4 qubits, 0(0) 
4 qubits, 0(8) 
5 qubits, 0(27) 


0.0589(62) 
0.0670(28) 
0.0267(89) 


0.01437(44) 
0.01325(16) 
0.01499(72) 


0.0352(31) 
0.0708(61) 
0.0254(23) 



(29) 



The statistical uncertainties of the fit values are given in parenthesis. They refer to the last 
two digits of the respective fit value. 
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Finally, we may again ask which kind of noise, static or random, is leads to a larger loss of 
fidelity. While in Fig. [3] (Larmor frequency noise) it was the random noise, here in Fig. Othe 
situation is more or less undecided. We note however that in the case of static noise (panal 
(a)) the different cases lead to a considerable spreading between the results in contrast to 
random noise (panel (b)), where the corresponding curves are very close together. We may 
also note that the 5-qubit case in panel (a) shows initially the fastest drop in fidelity. It 
therefore has the correct tendency if we believe that for more qubits and more Grover steps, 
the static errors should lead to lower fidelities [T9l [28] . 

Using the near resonant approximation, the numerical studies on Rabi frequency noise can 
be performed very efficiently. This would allow us to consider the Grover search algorithm in 
a larger quantum register with more qubits. We would need longer pulse sequences and more 
complex gate sequences for the implementation of the oracle. However, the main obstacle in 
considering more qubits are the unitary errors of our pulse sequence. Without additional error 
reduction {e.g. along the lines of Ref. [29]), the unitary errors alone would already reduce the 
final fidelity Fgnd so much that the effects of additional noise could not be analysed. 

6 Conclusion 

We studied the effect of noise on the performance of Grover's search algorithm implemented 
on a nuclear spin quantum computer with four spins coupled via first- and second-neighbor 
Ising interactions. Starting from the ground state, the quantum search algorithm attempts 
to build-up the target state in the quantum register, using the information obtained from 
inquieries of a quantum oracle. We used the fidelity to quantify the effect of noise on the 
final state, i.e. the result of the algorithm. We considered different types of noise: (i) Noise 
affecting the Larmor frequencies which would be due to fluctuations in the static magnetic 
field, and (ii) noise affecting the Rabi frequency, which would be due to fluctuations in the 
intensity of the radio-frequency pulses. We simulate the noise by randomly detuning the 
respective frequencies from their nominal values. For both cases (i) and (ii) we considered 
two types of noise, systematic and random errors. In the case of systematic errors, the 
resepective frequencies are randomly chosen at the beginning of the algorithm and then kept 
fixed during the whole protocol. By contrast, in the case of random errors, the respective 
frequencies are randomly chosen for each radiofrequency pulse. 

We found that the threshold of this amplitude for still having good performance of the 
algorithm is < 0.005 for Larmor-frequency noise and < 0.1 for systematic error on Rabi- 
frequency .We also found that this algorithm is much more sensitive to systematic errors than 
random errors which behave like an exponential and Gaussian decay respectively. One may say 
that this Gaussian behavior corresponds to the robustness of the algorithm. As we mentioned 
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before, computing time increases very much with the number of qubits, making the increasing 
of the number of Grover steps not feasible within our model. 
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